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Raman selection rules for electronic and magnetic excitations in BaFe2As2 were the- 
oretically investigated and applied them to the separate detection of the nodal and 
anti-nodal gap excitations at the spin density wave (SDW) transition and the separate 
detection of the nearest and the next nearest neighbor exchange interaction energies. 
Raman spectra are composed of magnetic excitations with gradually decreasing inten- 
sity toward far above the SDW transition temperature (Tgow) an d electronic excitations 
induced by the Brillouin zone folding below T SDW . The SDW gap has Dirac nodes, be- 
cause many orbitals participate in the electronic states near the Fermi energy. Using a 
two-orbital band model the electronic excitations near the Dirac node and the anti-node 
are found to have different symmetries. Applying the symmetry difference to Raman 
scattering the nodal and anti-nodal electronic excitations are separately obtained. The 
low-energy spectra from the anti-nodal region have critical fluctuation just above Tsdw 
and change into the gap structure by the first order transition at Tsdw> while those 
from the nodal region gradually change into the SDW state. Magnetic excitations are 
observed as a very broad peak in all polarization configurations. The selection rule for 
two-magnon scattering from the stripe spin structure was obtained. Applying it to the 
two-magnon Raman spectra it is found that the magnetic exchange interaction energies 
are not presented by the short-range superexchange model, but the second derivative 
of the total energy of the stripe spin structure with respect to the moment directions. 
The selection rule and the peak energy are expressed by the two-magnon scattering 
process in an insulator, but the large spectral weight above twice the maximum spin 
wave energy is difficult to explain by the decayed spin wave. It may be explained by 
the electronic scattering of itinerant carriers with the magnetic self-energy in the lo- 
calized spin picture or the particle-hole excitation model in the itinerant spin picture. 
The magnetic scattering spectra are compared to the insulating and metallic cuprate 
superconductors whose spins are believed to be localized. 
KEYWORDS: BaFe 2 As 2 , SDW gap, Dirac node, Two-magnon Raman scattering 

1. Introduction 

Two-dimensional iron pnictides form a new family of superconductors with the tran- 
sition temperature up to T c = 55 K. 1 ' 2 The BaFe 2 As 2 family undergoes superconducting 
state when the spin density wave (SDW) is suppressed by substituting the element 3-5 
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or pressurizing. 6 Quasi-static magnetic measurement shows the disappearance of the 
spin order and magnetic excitations above the spin density wave transition tempera- 
ture T S dw or m the metallic phase. However, the lost spin order is only perpendicular to 
the two-dimensional layer 7 and the high-energy magnetic excitations in the layer clearly 
remains above Tsdw or in the metallic phase as observed in magnetic Raman scattering 
and neutron scattering. 8-12 It strongly suggests that the superconductivity is induced 
by the same magnetic correlation as in the SDW state. 

The magnetism of iron pnictide cannot be simply understood by the localized or 
itinerant spin model, because the magnetic moment M = 0.8 ~ 0.9/xe/Fe is inter- 
mediate. 10, 13-18 The SDW phase is usually stabilized by the decrease of the electronic 
energy by opening of the gap at the Fermi energy in the folded bands. In the iron 
pnictide the SDW is not simply explained by the energy gain due to the gap formation 
induced by the Fermi surface nesting between the hole pocket at the T point and the 
electron pocket at the M point in k space. The spin structure is consistent with the 
Fermi surface nesting in BaFe2As2, 19-23 but not in FeTe, because the spin stripe direc- 
tion is different by 45 . 24 ' 25 Johannes and Mazin 16 pointed out that the energy gain is 
mainly one-electron energy balance for the antiferromagnetic (AF) spin patterns. The 
energy gain by the nesting is only a part of it. In fact a full SDW gap does not open 
in BaFe 2 As 2 as known from the fact that the electric conductivity in the SDW state is 
even better than in the normal state. It is caused by the remaining Dirac nodes. 26-29 
Furthermore the magnetic exchange interaction energies are not simply obtained from 
the conventional short-range superexchange mechanism. 30-34 The exchange interactions 
are very different, if they are obtained from the total energy of the long-range stripe 
spin structure. The exchange interaction in the y direction is antiferromagnetic in the 
former case, but ferromagnetic in the latter case. 

BaFe 2 As 2 undergoes the SDW state by the first order transition 9 ' 17 ' 35-39 accom- 
panied by the tetragonal (J4/mmm)-orthorhombic [Fmmrn) structural transition at 
Tsdw = 137 K. 4 ' 17 ' 40 The magnetic order in the SDW state is a stripe type in which 
nearest-neighbor spins are antiparallel in the x direction and parallel in the y direc- 
tion. 10 ' 13 ' 14 ' 30,40 ' 41 In this article we use the a and b axes as the tetragonal crystallo- 
graphy axes and the x and y axes as the orthorhombic axes which are rotated by 45° 
from the a and b axes. 

The electronic states in the SDW state is different from the conventional one. Usually 
the SDW state is insulating, because a full SDW gap opens at the Fermi energy. Despite 
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the SDW state, BaFe 2 As 2 is metallic because electronic states remain at the Fermi 
energy (E F ). It is expressed by the selective inter-orbital coupling in the multi-orbital 
bands of BaFe 2 As2. 20 ' 42-44 Wang et al. 26 and Ran et al. 27 predicted the nodal gaps 
composed of point-contact Dirac corns. The theoretical works of the Dirac nodes were 
reported 26-29 ' 45 ' 46 and some of them showed that the Dirac nodes are protected. 26-29 
Angle-resolved photoemission spectroscopy (ARPES) observed the Dirac nodes near 
E F . 47 The /c-linear dispersion of the Dirac corn is expressed by the massless relativistic 
Weyl equation 29 which expresses the motion of a neutrino. It is known that the two- 
dimensional massless Dirac Fermion exhibits a variety of unusual phenomena. 48 For 
example, the interaction to the acoustic phonon is very weak, causing a very large 
mobility. 

In order to find out whether the gap excitations near the Dirac node and the anti- 
node were separately detected, the orbital components in the electron and the hole 
bands were calculated in the two-dimensional two-orbital tight band model. The most 
area of the band in the k x -k y space is composed of the mixed xz and yz orbitals, but 
the electronic states along four lines are composed of pure xz or yz orbital. The same 
orbital components in two bands couple to open the SDW gap. The optical transition 
at the anti-node is from pure xz to pure xz or from pure yz to pure yz, while that at 
the Dirac node is from pure xz to pure yz or from pure yz to pure xz. The symmetry of 
the electronic transition at the anti-node is Ai g +B 2g and that at the node is B lg . The 
symmetry difference enables us to detect the anti-node and the node separately. The 
experimentally obtained B 2g spectra show critical fluctuation just above Tsdw and jump 
into the gap structure at Tsdw, while the B\ g spectra gradually change as temperature 
decreases through the Tsdw- 49 ' 50 The gap structure is composed of the reduced intensity 
below 300 cm -1 and the peaks at 400 and 800 cm -1 . The anti-nodal gap energies are 
the same as those in infrared spectroscopy 51 and ARPES. 52-57 

The magnetic exchange interaction energies also cannot be simply determined. The 
exchange interaction energies calculated in the superexchange interaction model are an- 
tiferromagnetic with the similar magnitude in the x and y directions, J\ x ~ J± y , because 
the anisotropy between x and y is small even in the orthorhombic SDW structure. 31 ' 33 
In order to stabilize the stripe spin structure, the next nearest neighbor exchange inter- 
action energy J 2 must be larger than half of the nearest neighbor exchange interaction 
energy. On the other hand the exchange interaction energies obtained from the second 
derivative of the total energy of the stripe spin structure with respect to the angle of 
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the moments are antiferromagnetic in the x direction and weak ferromagnetic in the 
y direction.. 32 ' 34 The superexchange interaction is related to the strong correlation of 
on-site Hubbard Coulomb repulsion. It was argued that the systems are moderately cor- 
related and the largely local iron moments are driven by Hund's intra-atomic exchange 
rather than by the on-site Hubbard repulsion. 16 

The magnetic excitations observed by neutron scattering are analyzed (1) by the 
antiferromagnetic interactions in both x and y directions 9 ' 14 ' 18 ' 30,41 ' 58 and (2) antifer- 
romagnetic in x and ferromagnetic in y. 14 ' 59 Usually the full spin wave dispersion is 
obtained in the localized spin picture. The similar dispersion is also obtained in the 
itinerant spin picture, but the high-energy parts are decayed into particle-hole excita- 
tions (Stoner continuum) and the continuum excitations extend above the maximum 
of the spin wave dispersion. 60-63 Many low-energy neutron scattering experiments did 
not reach the top of the dispersion. 9, 10 ' 14,30 ' 41,58 Diallo et al. 18 observed broad magnetic 
dispersion above 100 meV and suggested the damping of spin waves by the particle-hole 
excitations, while Zhao et al. 59 reported the full spin wave dispersion expressed by the 
antiferromagnetic in x and ferromagnetic in y. No consensus is obtained whether the 
itinerant or localized view is appropriate at present. 

In order to clarify whether the exchange interaction energies are determined by the 
short-range superexchange interaction mechanism or the total energy of the long-range 
stripe spin structure, the selection rule of two-magnon Raman scattering from the stripe 
spin structure was investigated. Two-magnon scattering in the insulating antiferromag- 
net is caused by the change of two spins AS Z = ±1. We found that the two-magnon 
scattering mechanism can be separately considered for the change of spins at the near- 
est neighbor sites and the next nearest neighbor sites. They give different symmetries 
and peak energies. Applying this selection rule to the two-magnon Raman scattering, 
it is found that the exchange interaction are given by the total energy mechanism. 

The peak energy of the present metallic magnetic Raman spectra is the same as the 
two-magnon scattering from localized spins in insulator, but the large spectral weight 
above twice the maximum energy of the spin wave dispersion cannot be understood 
by the two-magnon scattering process. Magnetic Raman spectra in metal is usually 
expressed by the two-magnon scattering from decayed spin waves. The decrease of the 
spin correlation length reduces the magnetic excitation energy as well as the life time. 
Therefore the large spectral weight at high energy cannot be interpreted by this model. 
We have to take into account the magnetic component of the itinerant carriers travelling 
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in the antiferromagnetic spin sea. The effect is interpreted both in the localized spin 
model and the itinerant spin model. In the localized spin model the high-energy mag- 
netic scattering is interpreted as electronic scattering of itinerant carriers in the AF spin 
lattice, because the moving carriers overturn spins and the electron spectral function 
have the magnetic component expressed by the self-energy in the string model. 64-71 In 
the itinerant model the magnetic scattering is interpreted as particle-hole excitations 
from the majority spin states to the minority spin states. 60-63 The magnetic scattering 
in metallic BaFe 2 As 2 is compared to the hole-doped La 2 _ x Sr :E Cu04 and electron-doped 
Nd2_ x Ce x Cu0 4 . The electron density dependence in BaFe 2 As 2 is found to be different 
from the cuprates 

The present experiments are based on the theoretical findings that the electronic 
excitations near the Dirac node and the anti-node can be separately detected using the 
different Raman symmetries and the nearest and the next nearest neighbor exchange 
interaction energies can be also separately detected using the different Raman symme- 
tries. Hence the relevant theories are presented followed by the experimental results. 
Section 2.1 presents the electronic Raman scattering mechanism of the SDW gap which 
is different from the normal metal and the superconducting gap. Section 2.2 presents 
the formulation of the nodal gap and anti-nodal gap in the two-orbital tight binding 
model. Section 2.3 presents the experimental procedure. Section 2.4 presents the exper- 
imental results of the high-energy spectra of the electronic Raman scattering and the 
magnetic Raman scattering. Section 2.5 presents the low-energy spectra of the Dirac 
node and the anti-node. Section 2.6 presents phonons. Section 3.1 presents the robust 
antiferromagnetic correlation in two-dimensional magnetism. Section 3.2 presents the 
symmetry of two-magnon scattering. Section 3.3 presents the experimental results of 
two-magnon Raman scattering. Section 3.4 presents the effect of conductive carriers to 
the magnetic scattering in comparison to the cuprate superconductors. 

2. SDW State 

2.1 Raman scattering of the SDW gap 

Sometimes the SDW gap is treated in the similar way as the superconducting gap in 
the Raman scattering experiment. However, their excitation processes are completely 
different. The Raman scattering of the superconducting gap is expressed by the intra- 
band transition with the momentum shift |fef — fei| = A A; at ki « kp, where kp is 
the Fermi wave vector. On the other hand the SDW gap scattering is expressed by the 
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Fig. 1. (color online) (a) Intraband electronic Raman excitations (red arrows) in the normal state. 
The green band has the electron pocket at k = ir and the blue band has the hole pocket at k = 0. 
The momentum shift in the Raman process is kept to be small in order to satisfy the momentum 
conservation with incident and scattered light, (b) In the superconducting state the excitation across 
the gap gives the pair breaking spectra. The electronic transition is limited to near the original Fermi 
surface. The initial and the final states are on the same band, (c) In the SDW state the electronic 
bands are folded into — tt/2 < k < tt/2. The intraband electronic excitations are limited to near the 
band crossing parts, while the interband transitions (purple dashed arrows) are available anywhere. 
Therefore the interband transition is dominant. 



interband transition between the original band and the folded band produced by the 
momentum shift of fe nest , where fe nest is the nesting vector (vr, 0). The vertical transition 
is allowed all over the folded Brillouin zone. 

Electronic Raman scattering in the normal state is induced by the first order of the 
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A 2 term and the second order of the p ■ A term in the electron-radiation interaction 
term (P — ~ C A) 2 . The matrix element is 72 ' 73 

1 ( sr^ < c > k i\Pp\b, h + q { X b, k, + q { \P a \a, k, > 

0a/3-\ > 



m = ( ;v 1 



1 s m 



+ x 



, (1) 



"7. \ ^ e t,fc+<Zi — e s,fc — 

where X is the term with the different time order, m the free electron mass, ef and ef 
polarization vectors of incident and scattered light, a and j3 the Cartesian coordinate, 
oui and q\ the incident photon energy and wave vector, a, 6, and c are the bands of the 
initial, intermediate, and final electronic states, and fcj and k { are the initial and final 
wave vectors of the electron. The scattered photon wave vector q s is fcj + q { — k { . The 
qi and q s are approximately zero. 

If the electron returns to the initial band "intraband transition, c = a" , eq. (1) is 
the same form as the k ■ p perturbation. In the low energy and the long wavelength 
approximation of the incident light eq. (1) becomes 

M = e?ed—\ , (2) 

where 1/m* is the effective inverse mass tensor at k^. 

If the final band c is different from the initial band a "interband transition, c ^ 
a", the scattering probability is determined by the symmetries of the initial and final 
electron wave functions. 

Figure 1 shows the electronic Raman excitations in (a) the normal metal, (b) super- 
conducting state, and (c) SDW state. In the normal state the momentum shift across 
the Fermi wave vector is the order of Ak < /czb/1000 to satisfy the momentum conser- 
vation with incident and scattered light, where kzB is the wave vector at the Brillouin 
zone boundary. The energy shift in Raman scattering from the intraband transition (red 
arrows) is usually limited to less than a few tens cm^ 1 , if the strong correlation effect is 
not taken into account. The scattering probability is expressed using the inverse mass 
tensor of eq. (2). 

The electronic transition in the superconducting state is the intraband transition. 
The electronic excitation is limited to near the gap because of the momentum conserva- 
tion with light. The scattering probability is expressed using the inverse mass tensor of 
eq. 

( 2 ) n, 74-76 The 

superconducting gap symmetry with the largest scattering intensity 
in BaFei.84Coo.i6As2 is f^g- 11 ' 50 ' 77 The scattering intensities from the electron pockets 
and hole pockets were calculated using the inverse mass tensors and the screening terms 
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in two-band model 11 and the full-potential linearized augmented plane wave (LAPW) 
model. 76 The importance of the multi-orbital was pointed 11,78 ' 79 

In the SDW state the electronic bands are folded into — it/ 2 < k < n/2 for the 
nesting vector fc nest = (vr,0). The intraband transition is limited to near the gap as in 
the superconducting gap. However, the dominant electronic excitations are the interband 
transition (purple dashed arrows) between the original band and the folded band all 
over the folded Brillouin zone. Among them the transitions at the rather wide area near 
the gap make the gap structure. The Raman intensity of the SDW gap is one hundred 
times larger than that of the superconducting gap. The electronic Raman scattering of 
the interband transition is not given by the inverse mass tensor. The symmetries of the 
wave functions directly result in the scattering probability. Then we have to know the 
nested band structure to understand the Raman scattering of the SDW gap. 

2.2 SDW state in the two-band model 

The electronic states near the Fermi energy is created from the two-dimensional FeAs 
layer. Wang et al. 2& and Ran et al. 27 pointed out that the SDW gap in BaFe 2 As 2 is not 
a full gap, but a nodal gap. In the two-orbital band model the hole pocket changes the 
dominant orbitals xz-yz-xz-yz on circling the Fermi surface, while the electron Fermi 
surface is represented by only one orbital. Under the conditions of (1) the collinear 
SDW order (denoted by the magnetization direction of x), (2) the inversion symmetry, 
and (3) the time-reversal symmetry, the vorticity is ±2 for the hole pockets and for 
the electron pockets. It creates the Dirac nodes with the vorticity ±1 in the SDW gap. 
The nodes are located at the intersection between the SDW wave vector and the hole 
Fermi surface in the two-band model. The Dirac nodes are shown to be stable in the 
five band model, too. 27 ' 29 

The different optical selection rule near the node and the anti-node is expressed 
by the two-orbital model of xz and yz. Using the transfer integrals by Qi et al. 42 and 
Raghu et al. 43 the tight-binding Hamiltonian is 




(3) 



where 2x2 matrix K{k) is 



K(k) = (e+(k) - h)t + e xy (k)ri + e_(fc)r 3 , 



e + (k) 



(ti + £2) (cos k x + cos k y ) — 4t 3 cos k x cos k. 
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€-(k) — — (ti — £2) (cos k x — cos k y ), 



€ 



■I'll 



(k) = — 4t 4 sin k x sin k y , (4) 



where d xza and d xz>a are the creation and annihilation operators of a d electron with 
the xz orbital and spin a, \i the chemical potential, are Pauli matrices, and r is the 
unit matrix. t\ is the transfer integral between xz orbitals and t 2 between yz orbitals 
at the nearest neighbor sites in the x direction. t 3 is the transfer integral between xz 
orbitals and t 4 between xz and yz orbitals at the diagonal sites. The band dispersions 
are given by 

E ± (k) = e + (k) ± v / e 2 (fc) + e ^(fc) - ix. (5) 

The upper and lower bands are called the electron band and the hole band, because 
they make the electron pockets and the hole pockets, respectively. The eigen functions 
are 



e xy{k) 



(6) 



where the normalization factors A± are 



ely(k) + el{k) ± ^(k)^el(k) + el y (k)] . 



(7) 



The xz and yz components on the electron and hole bands are given by eq. (6). 

Figure 2(a) shows the Fermi surfaces of hole pockets a and (5 and electron pockets 7 
and 5 in the Brillouin zone of one Fe atom per unit cell. The transfer integrals are ti = 
— 1, t 2 = 1.3, £3 = £4 = —0.85 in the unit of |£i|. 43 The gray arrows are unit vectors in the 
direction of (—e xy (k), e_(fe)), where e xy (k) represents the mixing between xz and yz 
orbitals and e_(fe) represents the difference of xz and yz component. 42-44 ' 85 The upward 
component of the gray arrow is the xz component (red) and the downward component 
is the yz component (green) for the electron band and opposite for the hole band. The 
hole Fermi surface is composed of two orbitals whose weights gradually change on going 
round the Fermi surface. The Fermi surfaces in the k x and k y directions from (0, 0) are 
composed of the pure xz and yz orbitals, respectively. The SDW state is formed by 



J. Phys. Soc. Jpn. 



FULL PAPERS 




Fig. 2. (color online) (a) Hole Fermi surfaces a and /? and electron Fermi surfaces 7 and 8. The xz 
and yz orbital components are shown by red and green according to the color scale. The gray arrows 
indicate the unit vectors parallel to (—e xy (k), e_(fe)). 42 - 43 The dashed lines shows the unit cell in the 
SDW state, (b) Dispersions of the electron and hole bands. The Ep is shown by the translucent blue 
plane, (c) Folded bands in the SDW state, (d) Nodal gapped bands near the Ep. The uppermost 
dispersion is removed, (e) The same as above in the expanded energy scale. Note that the front side 
is E — k x in (b), while E — k y in (c)-(e). (f)-(h) Dispersions in (k x ,Q), (k x ,ir), and (0, k y ). The circle 
denotes the Dirac node and the square denotes the anti-node. 



shifting the Brillouin zone by fc ncst = (vr,0). The reduced zone into — 7r/2 < k x < tt/2 
is shown by the dashed lines. Figure 2(b) shows the dispersions of the electron band 
and the hole band. The unit of the energy is \t\\. The weight of the xz and yz orbitals 
are represented by the scale from red to green. The translucent blue plane shows the 
E-p. Figure 2(c) shows the dispersion in a half of the folded Brillouin zone without the 
band mixing. The crossing parts may split by the mixing. Note that the front side of 
the cross section is the E — k v . 
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The local interaction to form the SDW with the wave vector (ir, 0) is 27 

#sdw = M ab J2(-l) ix (4,a^ iM -dl ai d iH ), (8) 

i 

where 

M ab = [(j) T + fan + (f) 2 T 2 + faTs]*, (9) 

and a and b are orbitals and i x is the x coordinate of the iron atomic site % in the unit 
of the inter-iron atomic distance. The (— l) lx term gives the twofold periodicity in the 
x direction. 0j is the mean field parameter to minimize the free energy including the 
correlation term. If the intra-orbital interaction r or r 3 is dominant, the Dirac node 
appears at k y = and 7r which are the same as the calculations by Ran et al., 27 Morinari 
et al., 29 and Kaneshita et a/.. 62,84 

Figures 2(d) and 2(e) show the gapped two bands near the E F . The same orbital 
components in two bands couple to open a gap, if the r or r 3 is dominant. The difference 
is scarcely found whether r or r 3 is dominant. The connecting point of two Dirac corns 
is a little below E F near (0, 0) and above E F near (0,7r), which causes small electron 
Fermi surfaces near (0, 0) and hole Fermi surfaces near (0, ir). The electronic transitions 
near the Dirac nodes appear at the energy near zero and those at the anti-node appear 
at 2Ai and 2A 2 . The gap energies are 2A X = 0.49 and 2A 2 = 0.39, when the coupling 
constant 0o — 0.3 or 3 = 0.3 in the energy unit of |ti|. The transitions near the Dirac 
node and anti-node have different symmetries, because the former is the transition 
between yz and xz while the latter between yz and yz near (0, 0) or xz and xz near 
(7r, 0). The Ak ~ transition e\ from yz to xz + yz at (0, 0) and €2 from xz to xz + yz 
at (0,7r) in Fig. 2(d) have large joint density of states. The energies are not strongly 
affected by the coupling term eq. (9). The dispersions in (k x , 0), (k x , n), and (0, k y ) are 
shown in Figs. 2(f)-(h). 



2. 3 Experimental Procedure 

Single crystals of BaFe 2 As 2 were grown by the self-flux method. Raman spectra were 
measured on the fresh cleaved surfaces in a quasi-back scattering configuration using 
5145 A laser light. The polarization configuration of the Raman spectra is presented by 
{e- x e s ). The crystallographic axes of the tetragonal structure are a and b. The bisecting 
directions are x and y. The x and y directions are rotated by 45° from the a and b. The 
observed Raman spectra include phononic, magnetic, and electronic excitations. Raman 
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active phonon modes are lA ig + lB lg + 2E g in the tetragonal structure. 80-83 The or- 
thorhombic structure (Fmmm) is the sub-group of the tetragonal structure (I4/mmm). 
The second order phase transition is allowed from the change of the symmetry, but 
whether the transition is the first order or the second order is determined by the free 
energy. It is known that the transition becomes the first order, if the SDW transition is 
accompanied by the structural phase transition. In the present case the first order tran- 
sition of BaFe2As2 has been confirmed in many experiments. 9 ' 17 ' 35-39 The (aa) spectra 
allow the A lg + B lg (A g + B lg ) modes, (ab) B 2g (A g ) mode, (xx) A lg + B 2g (A g ), and 
(xy) B lg (-Big) in the tetragonal (orthorhombic) structure. The Raman system was cal- 
ibrated using a black body radiation so that the intensity is proportional to the Raman 
probability independently to the scattered light wave number. 

2.4 High-energy spectra of the anti-nodal gaps and the transitions at (0, 0) and (0, 7r) 
points 

Figure 3 shows the temperature dependence of polarized Raman spectra. The sharp 
peaks at 181 and 215 cm -1 are caused by the A lg and B lg phonons, respectively. The 
broad hump with the top energy 2400 ~ 2500 cm -1 at 300 K is the two-magnon scat- 
tering peak. The two-magnon intensity increases with decreasing temperature. Below 
T S dw ~ 130 K new peaks appear at 400, 800, 2150, and 3500 cm -1 . 

As stated in § 2.1 the electronic scattering in the normal phase is limited to less 
than a few tens cm -1 . The wide-energy intensity is derived from only magnetic excita- 
tions. In order to estimate the interband transition activated by the SDW transition, 
the differential spectra between the SDW state at 10 K and the normal state at 200 
K are shown in Fig. 4(a). The intensity scales are the same for different polarization 
configurations. It is noted that the differential spectra include the electronic scattering 
component which is activated in the SDW state and the two-magnon scattering com- 
ponent whose intensity increases with decreasing temperature. In order to clearly show 
the magnetic scattering component, the scattering intensities in the normal states at 
200 K and 300 K are enlarged to the estimated level at 10 K in Fig. 4 (b). The zero 
levels are adjusted. The spectra at 200 K and 300 K are almost the same. It confirms 
that the spectral shape of the magnetic scattering is the same except the intensity in 
the normal phase. The extra humps appeared at 10 K are created by the interband 
electronic transitions. The differential spectra between 10 K and 200 K are shown by 
the black curve. The peak positions are shown by the arrows. The 800 and 400 cm -1 
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Fig. 3. (color online) Temperature dependence of polarized Raman spectra in BaFe2As2- 



peaks are assigned to the anti-nodal SDW gaps of 2Ai and 2A 2 and the 2150 and 3500 
cm -1 peaks to the transitions e\ at (0, 0) and e 2 at (0, ir) points as shown in Figs. 2(d), 
2(f), and 2(g). These peak intensities are large in (aa), (ab), and (xx) and small in (xy). 
The 2Ai and 2A 2 are the transition from yz to yz and from xz to xz, respectively. The 
ei is the transition from yz to yz and xz. The density of the final state is larger in yz 
than xz. The e 2 is the transition from xz to xz and yz. The density of the final state is 
larger in xz than yz. The symmetry of the transitions from xz to xz and from yz to yz 
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Fig. 4. (color online) (a) Differential Raman spectra between the SDW state at 10 K and the normal 
state at 200 K. (b) The Raman spectra at 10 K in the SDW phase and at 200 and 300 K in the normal 
phase. The 10 K spectra include electronic scattering activated by the Brillouin zone folding and 
magnetic scattering. The 200 and 300 K spectra is only magnetic scattering above a few tens cm -1 . 
The spectra at 200 K and 300 K are expanded in intensity and the zero levels are adjusted. The 
differential spectra between 10 K and 200 K are shown by the black curve, (c) Differential spectra 
obtained in the same way as (b). 



is A\ g + B 2g and that from xz to yz is B\ g . The (aa), (xx), (ab), and (xy) spectra are 
active in Ai g + Bi g , Ai g + B 2g , B 2g , and B\ g , respectively. Therefore the (xx), (ab), and 
(aa) spectra observe above four transitions strongly and (xy) spectra observe e± and e 2 
weakly. The experimental results roughly satisfy the selection rule. 

Figure 4(c) shows the differential spectra between the SDW phase and the normal 
phase at 200 K, and also between 200 K and 300 K in the normal phase. The 400, 800, 
2150, and 3500 cm -1 peaks start to increase below the Tsdw- The peaks are observed 



J. Phys. Soc. Jpn. 



FULL PAPERS 




at 150 K which is a little above Tsdw ~ 130 K of this sample. The intensities of these 
peaks increase like the second order transition as temperature decreases. 

Figure 5 shows the 10 K spectra at various polarization configurations and the 
optical conductivity at 10 K obtained by Hu et al.. 51 The 400 and 800 cm -1 SDW 
gap peaks in the present Raman spectra correspond to the 360, 890 cm -1 peaks in the 
optical conductivity. 51 They are also correspond to the gaps observed in ARPES. 47 The 
fine peaks are noticeable in Raman and infrared spectra. The 662 and 780 cm -1 peaks 
in the (xx) and (xy) spectra and the 882 and 1072 cm -1 peaks in the (aa) and (ab) 
spectra correspond to the fine structure in the optical conductivity spectra. 51 There 
are two possibilities for the origin of the fine structure. One is the electronic transition 
and the other is the phonon side band. Many other transitions may appear in the 
five band model. 27,84 Yi et al. 56 reported in ARPES that the electronic structures are 
reconstructed in the SDW state and may not be described in the simple folding scenario 
by the correlation effect, because the introduction of Hubbard U in the local-density 
approximation to fit the band dispersion changes the structure at Ep. The continuum 
between 400 and 800 cm -1 may have such effects. The Fermi surfaces have rather 
large three-dimensional characters 85 and it may induce some complicated structure. On 
the other hand the energy difference 118 cm -1 and 190 cm -1 between the main peak 
and the fine peak is close to the 117 cm -1 E g phonon energy and the 182 cm -1 Ai g 
phonon energy. 81 Those phonon modes forming the phonon side bands may be localized 
modes. 86 

The humps at 2150 and 3500 cm" 1 are not observed in the infrared spectra, because 
the d-d transition is infrared inactive. We assigned the humps to the ei and transitions 



J. Phys. Soc. Jpn. 



FULL PAPERS 




Fig. 6. (color online) Temperature dependent low-energy Raman spectra in BaFe2As2- The low- 
energy spectra below 300 cm -1 jump to the high-intensity level at Tsdw and then gradually decreases 
as temperature increases in the (xx) and (ab) spectra. The red arrows show the jumps. The (eta) and 
(xy) spectra do not show the jump. 



in Fig. 2(d). Those transitions may be related to the flat bands whose energies decrease 
with the increase of U. The U reduces the magnetic moment M = 0.5 ~ 1 /xb/Fg in the 
local-density approximation plus Hubbard U (LDA+U) band calculation. 56 However, it 
is difficult to compare the Raman peaks to the band calculation, because the bandwidth 
is renormalized and the energy shift is different at the F and X points to fit the ARPES 
data. 

2.5 Low-energy spectra of the Dirac node and the anti-node 

The temperature dependence of the polarized low-energy Raman spectra is shown 
in Fig. 6. All the spectra are plotted in the same intensity scale. The sharp peaks at 
181 and 215 cm -1 are the A\ z and B\ g phonons in the tetragonal phase, respectively. 
They are discussed later. The 400 and 800 cm -1 humps are observed at 10 K in all 
polarization configurations. However, the temperature dependence below 300 cm -1 is 
very different between the (aa) and (xy) group and the (xx) and (ab) group. In the 
(xx) and (ab) spectra the low-energy scattering intensity abruptly increases from 125 
K to 130 K in the (xx) spectra and from 120 K to 125 K in the (ab) spectra, and then 
gradually decreases as temperature increases. 49,50 On the other hand the jump is not 
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Fig. 7. (color online) Differential spectra from 200 K. The intensity below 200 cm" 1 increases in 
the (xx) and (ab) spectra as temperature decreases below 150 K and jumps to the low intensity at 
Isdw- The SDW- induced anti-nodal gap peaks at 400 and 800 cm" 1 gradually increase in intensity as 
temperature decreases below a little above Tsdw m a U polarization configurations, but is weakest in 
(xy). The zero levels are shown by the same color lines. 

observed in the (aa) and (xy) spectra near the Tg DW . The small temperature difference 
of the jump in (xx) and (ab) is induced by the first order transition. 

It is noted that the experiment was made in the sequence of (aa) and (ab) at 10 K, 
(aa) and (ab) at 100K, • • • , 300 K and after rotating the sample by 45° (xx) and (xy) 
at 10 K, (xx) and (xy) at 50 K, ■ ■ • , 300 K. The (aa) and (ab) spectra at the same 
temperature are obtained exactly in the same condition except for the scattered light 
polarization. The same is for (xx) and (xy). Therefore the spectral jump at Tgpjvv in 
the (xx) and (ab) spectra and the continuous change in the (aa) and (xy) spectra are 
intrinsic properties of this material. 

Figure 7 shows the differential Raman spectra between various temperatures and 
200 K. The zero levels are shown by the horizontal lines with the same color as the 
spectra. In the (xx) and (ab) spectra the intensity below about 200 cm -1 increases with 
decreasing temperature from 150 K to Tsdw ~ 130 K and then jumps to decrease by 
the opening of the SDW gap. On the other hand the jump is not observed in the (aa) 
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and (xy) spectra. 

The above difference between the (xx) and (ab) group and the (aa) and (xy) group 
is caused by the different symmetries for the electronic transitions near the anti-node 
and the Dirac node. The electronic transition near the anti-node is from xz to xz or 
from yz to yz and that near the Dirac node is from xz to yz or from yz to xz. The 
symmetry of the transition near the anti-node is A lg + B 2g and that near the node is 
.Big. The (xx) spectra are active to A lg + B 2g and (ab) to B 2g , (aa) to A lg + B lg , and 
the (xy) to B Xg . Therefore the (xx) and (ab) spectra represent the electronic transitions 
near the anti-node and the (aa) and (xy) spectra represent the Dirac node, if the A lg 
component from the anti-node is ignored. The spectra representing the anti-node have 
the jump into the gap spectra. On the other hand the spectra representing the Dirac 
node has not the specific change at T SDW - The A Xg component might be observed in 
the (aa) spectra, but the spectra do not show the jump at Tsdw- Instead the (aa) 
spectra gradually change into the spectra of the gapped excitations. In contrast to the 
(aa) spectra the gap spectra are weak in the (xy) spectra which represents only the 
excitation near the node. 

The increase of the low-energy scattering intensity on approaching the Tsdw from 
high temperature is induced by the critical fluctuation relating to the opening of the 
anti-nodal gap, because it is observed only in the B 2g symmetry in which the anti-nodal 
gap excitations are observed. The magnetic fluctuation above Tsdw was also observed 
in neutron scattering. 10 The critical fluctuation is the characteristic properties of the 
second order phase transition. The simultaneous SDW and structural transitions in 
BaFe2As2 splits into two second-order-like transitions on substituting Co for Fe. 36 The 
phase locking of the SDW to the lattice potential by the simultaneous structural tran- 
sition changes the transition into the first order and the SDW into the commensurate 
one. The SDW and structural phase transition of BaFe 2 As2 is, however, close to the 
second order transition. It induces the critical fluctuation. 

The above analysis is done on the assumption that the dominant orbital mixing 
term is r or r 3 in eq. (9). If T\ or r 2 is dominant, the positions of the Dirac node and 
the anti-node are exchanged and the observed spectral symmetries at the Tsdw become 
inconsistent to those of the electronic transitions. 

As discussed above the symmetry dependence of the observed spectra is consistently 
interpreted by the electronic excitations near the anti-node and the Dirac node. How- 
ever, the dominant orbitals at the nesting parts in electron pocket may be modified in 
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the five band model. 20 ' 44, 45 ' 85 ' 87 ' 88 More detailed theoretical analysis is expected. 

2.6 Phonon mode activated in the SDW state 

The sharp peaks in the tetragonal phase above T S dw are observed at 181 (Ai g ) 
cm -1 in the (aa) and (xx) spectra and 215 cm -1 (Bi g ) in the (aa) and (xy) spectra 
of Fig. 6. 80-83 Below T SDW the (A lg ) peak appears at 187 cm -1 in the (ab) spectra. 
The Ai g mode in the tetragonal (IA/mmm) structure changes into the A g mode in 
the orthorhombic structure by the compatibility relation of the group theory. Taking 
into account the 45° rotation of the crystallographic axes at the phase transition from 
the tetragonal to the orthorhombic structure, the Raman intensity of the A g mode 
in the (ab) (in the tetragonal notation) polarization configuration is proportional to 
\Ru — /?22 1 2 using the Raman tensor in the orthorhombic structure. The intensities in 
the (aa), (xx), and (yy) spectra are proportional to \\Ru + -R22I 2 , I-R11J 2 , and I-R22I 2 , 
respectively. The (xx) and (yy) spectra may mix if the crystal is twinned, but both 
spectra represents the same A g mode and the difference is small even in the untwinned 
crystal. Other polarization configurations are not affected by the twin structure. The 
intensities in the (aa), (xx), and (yy) spectra are expected to be nearly the same and 
much larger than that in the (ab) spectra because Ru ~ i? 22 . Instead at 10 K the 
intensities in the (xx) and (ab) spectra are 1.5 ~ 2 times as large as the intensity in 
(aa). This phonon is the mode in which As atom moves in the c direction. 80 It has large 
magneto-phonon interaction, because the Fe-As-Fe angle is very sensitive to the Fe-Fe 
exchange interaction energy. 31,32 ' 89 The detailed mechanism of the enhancement is still 
an open question. The enhancement of the A ig phonon below Tsdw was reported in 
CaFe 2 As 2 82 and BaFe 2 _ :r Co :r As 2 . 83 Similar enhancement of the infrared active phonon 
was reported in BaFe 2 As2 90 

3. Magnetic Raman Scattering 

3.1 Antiferromagnetic correlation in low- dimensional magnetism 

The exchange interaction energies are very different whether the calculation starts 
from the short-range superexchange interaction or from the total energy of long-range 
AF stripe spin structure. In the former case the calculation of the local superexchange 
interaction gives AF exchange interaction energies for all directions reflecting the equiv- 
alent x and y directions. The stripe spin order is stable if J\ x = J\ y < 2J 2 , otherwise 
the checkerboard spin order is stable. 31,33 On the other hand in the latter case the 
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exchange interaction energies are obtained from Jij(R) = —d 2 E/d9i(0)d9j(R), where 
E is the total energy and 9j(R) is the angle of the moment of the jih spin. 32,34 The 
calculated energies are antiferro magnetic SJi x = 43, ferromagnetic SJ ly = —3.1, and 
antiferromagnetic diagonal SJ 2 = 14.3 meV. 34 The low-energy spin wave in the SDW 
state was observed by neutron scattering in BaFe2As2, 10,14 SrFe2As2 30 and CaFe2As2- 41 
Ewings et al. 14 fitted the spin wave velocity at 7 K in BaFe2As2 by (1) the long-range 
exchange interaction model with SJ lx = 36, SJ ly = —7.2, and SJ 2 = 18 meV and 
(2) the short-range exchange interaction model with SJ lx = 17, 5, SJ ly = 17.5, and 
SJ 2 = 35 meV. Zhao et al. disclosed in high-energy neutron scattering that the entire 
spin wave dispersion of CaFe 2 As 2 (T N w 170 K) at 10 K is expressed by the long-range 
model with SJ lx = 49.9, SJ ly = -5.7, SJ 2 = 18.9, and SJ Z = 5.3 meV. 59 
The spin wave dispersion ignoring the anisotropy terms is given by 14 

hu{k) = VA 2 - D 2 , (10) 
+ J lx + 2J 2 + J z ), (11) 



A = 2S { J ly 



cos — — 1 
2 



D — 2S I J\ x cos — + 2 J 2 cos — cos — + J z cos k z 1 . (12) 



14 



Figure 8 shows the dispersion at k z = with (a) SJ lx = 36, SJ\ y = —7.2, SJ 2 = 18, 
and SJ Z = 5.3 meV, (b) SJ\ X = 17.5, SJ\ y = 17.5, SJ 2 = 35, 14 and SJ Z = 5.3 meV, 
and (c) SJi x = 40, SJi y = 40, SJ 2 = 20, and SJ Z = 5.3 meV. (a) and (b) have nearly 
the same spin wave velocity, but the maximum energy in (b) is lower than in (a). The 
energy at k x = or k y = goes to zero at the critical point of (c) J lx = Ji y = 2J 2 in 
Fig. 8(c). 

The superconducting phase does not overlap the SDW phase or slightly overlaps 
at the crossover region, if it is observed by quasi-static experiments such as magnetic 
susceptibility, resistivity, specific heat, /iSR, and crystal structure. 36, 91-93 Sometimes it is 
called as the evidence of competing superconductivity and magnetism. However, such a 
phase diagram is valid only at k and 00 ~ 0. Neutron scattering experiments in high 
temperature superconductors disclosed that only the magnetic correlation perpendicular 
to the CuC>2 plane disappears above the magnetic transition temperature. 8-10 This is 
a common property in low dimensional magnetism. The AF transition temperature Tn 
in two-dimensional antiferromagnet is determined by the weakest exchange interaction 
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J z cLS 

where L is the intralayer AF spin correlation length and d is the inter-spin length in 
the layer. Neutron scattering showed that the correlation length is 15 ~ 18 A at 136 
K (~ Tsdw) an d it remains to be 6 ~ 8 A at 180 K. 9,10 The magnetic excitations 
with the wavelength longer than the correlation length may disappear but those with 
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Fig. 9. (color online) Spin structures in (a) the checkerboard type in the cuprate superconductors 
and (b) the stripe type in BaFe2As2. The AF exchange interactions are presented by the line segments 
and the ferromagnetic interaction by the dashed line segment. Note that a and b are the tetragonal 
crystallographic axes at high temperatures. 



shorter wavelength survive above Tn and even in the superconducting phase. Two- 
magnon Raman scattering in the AF insulator represents the joint density of states of 
k and —k magnons to satisfy the momentum conservation with light. The dominant 
part is derived from the zone boundary with the large density of states. Therefore the 
two-magnon peak survives till the last as the correlation length decreases. 

3.2 Symmetry of two-magnon scattering induced by Coulomb interaction in insulator 
Two-magnon scattering in AF insulator arises from the simultaneous changes of 
AS Z = ±1 on both spin sublattices by the Coulomb interaction. 94-99 The two-magnon 
Hamiltonian can be written as 95 

-f^two-mag = ■^a^ys{ r ) e i' e B ^R^R+rJ (14) 

Rr 

a/?7<5 

where e- x and e s are incident and scattered light polarizations and r is the vector from 
the site R to the antiparallel spin site. B is determined so that H two ^ mag is totally 
symmetric. The spin term *Sjj*Sr, +1 . is taken so that the total z component of excited 
spins is zero as 

qX qX , qV qV _ 1 I C+ C~ I 0+ \ ( 1 Z\ 

•JR^R+r "T D i? D fi+r — ^ R R + r ^R^R+r)' \ L0 ) 

In order to clarify the difference in the selection rule of the stripe AF spin structure 
in BaFe2As2 and the checkerboard AF spin structure in cuprate superconductors, the 
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checkerboard-type is first discussed. The spin structure is shown in Fig. 9(a). The 
nearest neighbor spins are in the a and b directions. All the nearest neighbor spins 
are antiparallel. The next nearest neighbor spins are parallel and cannot contribute to 
two-magnon scattering. The next nearest neighbor exchange interaction is much smaller 
than the nearest neighbor exchange interaction. The system Hamiltonian is 

H = jJ2Sk-S l . (16) 

kl 

The two-magnon scattering Hamiltonian is given by 97 ~" 

#two-ma g = ^2 A ( e i ■ r k i){e s ■ r kl )(S k ■ Si), (17) 

kl 

where r k i is the unit vector connecting sites k and / sites. Here the spin longitudinal 
term S R S R+r is added to eq. (15). Two-magnon scattering is active in (aa) and (xy) 
and inactive in (ab). In (xx) the two-magnon scattering Hamiltonian is the same as 
the system Hamiltonian except for the proportionality constant and the two-magnon 
scattering is inactive, because the two-magnon Hamiltonian commutes to the system 
Hamiltonian. The two-magnon peak energy is lower than the peak energy in the density 
of independent two magnon states by the magnon-magnon interaction energy J, because 
the exchange interaction energy between the nearest neighbor overturned spins is the 
same as before. 

Figure 9(b) shows the stripe spin structure in BaFe 2 As 2 . It should be noted that 
the nearest neighbor directions are x and y in BaFe2As2 instead of a and b in cuprate. 
The a and b axes are the tetragonal crystallographic axes. The nearest neighbor spin ar- 
rangement has two-fold rotational symmetry, antiparallel in the x direction and parallel 
in the y direction. All the second nearest neighbor spins are antiparallel. The system 
Hamiltonian is 

H = Jl x ^ Sr ■ SR +r + Jly ^ Sr ■ SR +r 

R,r=±x R,r=±y 

+ -h S * S R+r (18) 

R,r=±x±y 

The nearest and the nest nearest neighbor spins should be taken into account, because 
both exchange interaction energies are the same order, for example, S J\ x = 36 meV and 
SJ 2 = 18 meV. 14 The two-magnon scattering Hamiltonian from the nearest neighbor 
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spins is given by 

-^two-mag = Cefe*(SxS R+r + S R S R+r ), (19) 

r=±x 

because only antiferromagnetic spins in the x direction can contribute to the scattering. 
Two-magnon scattering is active, if both incident and scattered light has the electric 
field component in the x direction. Then allowed polarizations are (xx), (aa), and (ab). 
The crystals have magnetic twin structure for the interchange of x and y axes. The 
(yy) and (xx) polarizations cannot be distinguished. Then two-magnon scattering is 
inactive only in (xy). The magnon-magnon interaction energy for the nearest neighbor 
spin exchange is J\ x . 

The two-magnon scattering Hamiltonian from the second nearest neighbor spins are 
the same as eq. (17) of the checkerboard spin structure. The two-magnon scattering is 
active in (xy) and (aa), and inactive in (ab). The (xx) polarization is active, because 
the two-magnon Hamiltonian no longer commutes to the system Hamiltonian. The 
magnon-magnon interaction energy for the second neighbor spin exchange is J 2 - 

The sum of the nearest and next nearest neighbor terms gives two-magnon inten- 
sity. Therefore it is active in all polarization configurations. In the case of antifer- 
romagnetic Ji x and ferromagnetic J\ y the two-magnon scattering energy without the 
magnon-magnon interaction is the same in all polarization configurations, because the 
spin wave energy keeps constant along k y = 2n in Fig. 8(a). The magnon-magnon inter- 
action reduces the peak energy from twice the maximum energy of the density of states 
by Ji x in (ab) and by J 2 in (xy), because only the nearest neighbor gives scattering in 
(ab) and the next nearest neighbor gives scattering in (xy). In the (aa) and (xx) spectra 
both components contribute. 

3.3 Two-magnon scattering in BaFe 2 As 2 

Figure 10 shows the polarized Raman spectra in the normal state at 300 K. The 
broad peak at 2500 cm^ 1 is the two-magnon peak. In order to compare the two-magnon 
peak energies at 300 K in different polarization configurations, the intensities are ad- 
justed so that the peak heights are the same. The spin wave densities of states cal- 
culated using eq. (10) with J\ x = 36, J\ y = —7.2, J 2 = 18, and J z = meV (black) 
and Ji x = 17.5, J\ y = 17.5, J 2 = 35, and J z = meV (blue) are shown in Fig. 10. 
The energy is doubled from that of the density of single spin wave states. It is the 
density of independent two magnon states. The peak energy should be higher than the 
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Fig. 10. (color online) Polarization configuration dependence of the two-magnon peak. The sharp 
curves are the density of independent two magnon states calculated using eq. (10) with the parameters 
J\ x = 36, Ji y = —7.2, J 2 = 18, and J z — meV (black) and J\ x = 17.5, J\ y = 17.5, Ji — 35, and 
J z = meV (blue). Two-magnon scattering peak energy is lower than the peak energy in the density 
of independent two magnon states by the exchange interaction energy between overturned spins. The 
(ab) and (xy) peak energies are smaller by J\ x and J2, respectively. The (aa) and (xx) are intermediate 
between them. 



two-magnon scattering peak energy, because the magnon-magnon interaction reduces 
the two-magnon scattering peak energy. Therefore the blue density of states which is 
obtained from the superexchange interaction model is ruled out. The peak energy in 
the calculated density of independent two magnon states (black) is 2790cm -1 . 

The experimentally obtained two-magnon peak energy in the (ab) spectra is lower 
than the independent two magnon peak energy (black) by Ji x and that in the (xy) 
spectra is lower by J 2- The reduced energy J\ x and J2 are the expected magnon-magnon 
interaction energies for the (ab) and (xy) spectra assuming S = 1 as discussed at the 
end of § 3.2. In the (aa) and (xx) spectra the decreased energy is the average of J\ x 
and J 2- The experimentally obtained peak energy of two-magnon scattering is perfectly 
expressed by the total energy model of the long-range stripe spin structure. The very 
broad peak may be attributed to the short correlation length in the metallic state. 
However, the large weight at the higher energy than twice the highest spin wave energy 
may not be explained by the above Coulomb interaction-induced model. It is discussed 
in the next subsection. 

Very recently Chen at al. calculated two-magnon scattering peak energies. 100 The 
Bi g peak energy is a little lower than the Ai g + i?2 g peak energies in the model with 
Ji y = —0.1Ji x and J 2 = Jix in agreement with our result. They also showed that the 
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Fig. 11. (color online) (a) Temperature dependent polarized Raman spectra in the paramagnetic 
phase of BaFe2As2 (upper panel) and the superconductor BaFei.84Coo.i6As2 (lower panel). The density 
of independent two magnon states is shown by the dashed line. The two-magnon peak energy and 
the independent two magnon excitation energy are shown by the dashed lines. The inset is the spin 
structure, (b) Two-magnon Raman spectra in the insulating La2CuC>4 (upper panel) and the optimally 
hole-doped (x = 0.15) La2- a; Sr a ;Cu04 (lower panel). The density of independent two magnon states 
using single J* — 118 meV which is obtained from the two-magnon Raman scattering is shown by the 
solid line and that using J* — 104 and J% = —18 meV which are obtained from neutron scattering 103 
is shown by the dashed line in the upper panel. The low-energy peak at x = 0.15 in the (ab) spectra 
of La2_ a; Sr 2 ;Cu04 is the superconducting pair breaking peak, (c) Magnetic Raman spectra in the 
insulating Nd2CuC>4 (upper panel) and the electron-doped superconducting Nd2- a; Ce. E Cu04 (lower 
panel). The density of independent two magnon states with J* — 855 cm -1 is shown by the solid line. 



Bi g peak energy with J2 = J\ x = J\ y is very low. The B\ g spectra observe the spin wave 
in k space with the weight of (cos k x — cos k y ) 2 . The energy decreases as J2 approaches 
Ji x /2 = J\ y /2 in Fig. 8(c). However, the Bi g peak energy may not be so small, because 
the neutron scattering experiment shows J2 ~ 2Ji x = 2Ji y in the case of J\ x = Jiy 14 

3.4 Magnetic excitations induced by moving carriers in the localized model 

The magnetic Raman scattering in cuprate superconductors are considered by 
the localized model. Figure 11 shows the comparison among (a) iron pnictides 
BaFe2- x Co a; As2, (b) hole-doped superconductor L^-xSr^CuC^, and (c) electron- 
doped superconductor Nd2_ x .Ce x Cu04. Figure 11(a) shows the antiferromagnetic metal 
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BaFe 2 As 2 with T S dw = 130 K and the superconductor BaFei. 8 4Coo.i6As 2 with T c = 25 
K. The sharp dashed peak in the upper panel is the calculated density of independent 
two magnon states with SJi x = 36, SJ ly = —7.2, and SJ 2 = 18 meV. The two-magnon 
peak remains even in electron-doped superconductor BaFe1.84Coo.i6-A.S2, indicating the 
short-range AF magnetic correlation remains. 11 The two-magnon peak is very broad 
and has large spectral weight above twice the maximum energy of the spin wave (right 
dashed line). The high-energy weight cannot be interpreted by the decayed magnon in 
the Coulomb interaction-induced two-magnon scattering process, because the higher- 
order multi-magnon scattering intensity is known to be weak. 101 

In La 2 Cu0 4 of Fig. 11(b) the sharp peaks from 700 to 1500 cm -1 are two-phonon 
scattering peaks. The two-magnon scattering induced by the Coulomb interaction pro- 
cess is active in the B Xg symmetry which is observed in the (xy) spectra. The rather 
sharp peak around 3206 cm^ 1 at 5 K in the insulating phase is the two-magnon peak. 
The B 2g scattering intensity in the (ab) spectra is weak. If only the nearest neigh- 
bor exchange interaction is taken into account, the two-magnon peak energy is 2.76J 
(3. 38 J*), 102 where J* is the exchange interaction energy with the correction of the quan- 
tum spin effect of S = 1/2 by Oguchi. The obtained J* is 949 cm" 1 (118 meV). The 
maximum energy in the density of independent two magnon states is 4 J (4. 63 J*). 102 The 
energies of the experimentally obtained two-magnon Raman peak and the independent 
two magnon peak are shown by the vertical dashed lines. The density of independent 
two magnon states calculated using single J* is shown by the solid line in the upper 
panel of Fig. 11(b). Coldea et a/. 103 obtained the full spin wave dispersion in neutron 
scattering. The fitted parameters are J* = 104 and J| = —18 meV, where J 2 is the 
second nearest neighbor exchange interaction energy. The density of independent two 
magnon states is shown by the dashed line. The two-magnon scattering peak and the 
independent two magnon peak widely separates in the cuprates compared to BaFe2As 2 , 
because the spins in cuprates are S — 1/2 and those in iron pnictide are S = 1. There is 
small spectral weight above the independent two magnon peak energy in the insulating 
phase, indicating that the higher order multi-magnon scattering intensity is small. 101 
The two-magnon peak remains as temperature increases above the AF transition tem- 
perature Tn. In hole-doped La 2 _ x Sr :c Cu04 the peak energy shifts to low energy with 
broadening. The magnetic scattering cannot be derived from the same process as in the 
insulating phase, because the integrated intensity does not decrease. If the decrease of 
the peak energy is caused by the shortened correlation length, the integrated intensity 
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decreases. 104 

In Nd 2 Cu0 4 of Fig. 11(c) the two-magnon spectra in the insulating phase are nearly 
the same as La 2 Cu0 4 . The peak energy is 2890 cm -1 and the obtained J* is 855 cm -1 
The density of independent two magnon states using only J* is shown by the solid curve 
in the upper panel in Fig. 11(c). The electron doping causes the abrupt shift of the peak 
energy to high energy at the insulator-metal transition. The peak energy little shifts as 
the electron density increases in the metallic phase. The high-energy shift in the metallic 
phase cannot be interpreted by the Coulomb interaction induced two-magnon scattering 
mechanism. The spectral weights do not decrease in the metallic phase similarly to 
La 2 _ a: Sr :E Cu04. The B 2g spectral intensity observed in the (ab) spectra increases to 
nearly the same magnitude of the Bi g intensity in the (xy) spectra as carrier density 
increases. The magnetic scattering spectra in the metallic phase may be interpreted by 
electronic scattering with the self-energy part of the magnetic component as discussed 
in the following. 

The electronic state of a hole moving in the AF lattice is extensively investigated 
in cuprates. Hopping of a hole is caused by the back-hopping of an electron with the 
opposite spin from the nearest neighbor site. Therefore the moving hole leaves behind an 
overturned spin trace in the AF spin lattice. The increased magnetic energy prevents 
the hole hopping. A pair of opposite spins on the trace may return to the original 
directions by the quantum spin fluctuation, one and a half rotation of a hole around 
the nearest-neighbor four Cu square sites, 105 the short range AF correlation, and the 
hopping of other holes. The electronic state was calculated, for example, by the string 
model taking into account the quantum spin fluctuation. 64-71 The hole sees a linearly 
rising potential on increasing the migration distance in the limit of t » J, which 
makes a series of discrete energy levels 

En /t = -b + a n (J/tfl\ (20) 

where a n are the eigen values of a dimensionless Airy equation, a\ = 2.16, a 2 = 5.46, 
a3 = 7.81, and b = 3.28 at (ir/2, 7r/2). 69 The discrete peak height decreases with in- 
creasing n and connects to a continuum at high energies. The ground-state E\ is a 
delta-function and the width of the second peak is about 2 J at J = O.lt The energy 
difference between E 1 and E 2 is 4700 cm -1 using the widely accepted values of t = 0.4 
eV and J/t = 0.3. 106 This energy is close to the broad peak energy in Nd 2 _ x Ce x Cu0 4 
with x = 0.14. The discrete levels cannot be detected in Raman scattering, because the 
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levels have dispersions. The nearly constant peak energy as the carrier density increases 
in the metallic phase of Nd 2 _. r Ce, E Cu04 is successfully explained by the string model. 

In the metallic phase of the hole doped cuprate La 2 _ x Sr,,.Cu04, the magnetic peak 
energy continuously decreases with increasing the carrier density. The mean free path 
in the underdoped phase of the hole doped cuprate is one order shorter than that of 
electron doped cuprate. The migration distance is not as long as to allow the linearly 
rising potential approximation. The AF correlation length decreases with the increase of 
the hole density. 107 The electron Green function has a self-energy part of excited spins 
by the hole hopping. The electron spectral functions presented by the imaginary part 
of the retarded Green function has the coherent component near the original electronic 
energy and the incoherent component near the spin excitation energies. The Ak ~ 
transition between the electron spectral functions gives the spectral weight at the spin 
excitation energies. It is in contrast to the uncorrelated electron spectral function which 
is presented by the delta function at the energy of the dispersion and the Ak ~ 
transition gives only Ae ~ excitations. 

The very broad two-magnon peak in BaFe 2 As 2 is compared to the magnetic exci- 
tation peak in the metallic cuprate. The hole density is 0.081 and the electron density 
is 0.069 per iron atom at 140 K in BaFe 2 As 2 . 56 The spectra little changes by elec- 
tron doping in BaFe 2 _ :c Co :E As 2 . The large spectral weight above the maximum of the 
independent two-magnon energy and the nearly independent peak energy on the elec- 
tron density are similar to Nd 2 _ :r Ce x Cu04. However, the peak energy is given by the 
two-magnon excitation energy estimated from the spin wave dispersion in BaFe 2 As 2 , in 
contrast to the much higher energy in Nd 2 _ x Ce :r Cu04. 

3.5 Particle-hole excitations in the itinerant model 

In the itinerant spin model the particle-hole continuum is created above the spin 
wave dispersion. 60-63 The spin wave dispersion strongly decays when the energy overlaps 
with the particle-hole continuum. Kaneshita and Tohyama 62 calculated the imaginary 
part of the susceptibility using the five-band model. The spin wave dispersion becomes 
clear as the magnetic moment increases. The full spin wave dispersion with the max- 
imum energy 0.25 eV and the dispersion with the maximum energy 0.4 eV above it 
are obtained at the magnetic moment M = 0.8/zs/Fe. The calculation showed that 
the integrated intensity over k increases from energy zero to 0.2 eV. The intensity is 
larger in the paramagnetic phase than in the AF phase. In the Raman scattering two 
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particle-two hole excitations are observed. The scattering intensity decreases on increas- 
ing temperature. The increase of intensity in the paramagnetic phase is not observed. 
The calculation did not show whether the intensity maximum occurs at twice the top 
of the spin wave dispersion or not. The imaginary part of the susceptibility for the lon- 
gitudinal spin excitations and the charge excitations has large intensity from 0.2 to 0.4 
eV. The two particle-two hole excitations may observe from 0.4 to 0.8 eV. The observed 
spectra have the broad peak as shown in Fig. 11(a). 

4. Conclusion 

The change of the electronic and magnetic states at the SDW transition in BaFe2As2 
was investigated by Raman scattering. The tight binding band model with two orbitals 
disclosed that the electronic excitations near the Dirac node and the anti-node in the 
SDW state have different symmetries. Utilizing this selection rule the electronic excita- 
tions near the Dirac node and the anti-node were separately obtained. The electronic 
excitations near the anti-node have critical fluctuation above T SDW and change into the 
gap structure by the first order transition, while those near the Dirac node gradually 
changes. As for the magnetic scattering the selection rule of two-magnon scattering in 
the stripe spin structure was obtained. The magnetic exchange interaction energies are 
not presented by the superexchange interaction model, but the second derivative of the 
total energy of the long-range stripe spin structure with respect to the moment direc- 
tions. The magnetic excitation peak survives far above T SDW and in the superconduc- 
tor phase of BaFe 2 - 2: Co x As 2 , indicating that the short-range spin correlation remains 
in the two-dimensional magnet. The two-magnon peak energy is well understood by 
the conventional two-magnon scattering process in insulator. However, the high-energy 
spectral weight above twice the maximum energy of the spin wave may be explained 
by the string model in the localized spin picture or the particle-hole excitations in the 
itinerant picture. The difference from the high temperature superconducting cuprates 
were discussed. 
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